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ABSTRACT 


This report presents the results of an experimental 
investigation of the behavior of the turbulent early wake 
behind three flat plates of varying trailing edge config- 
uration submerged in water. 


The circumstances surrounding "singing" propellers and 
the influence of vortex formation upon this and other 
hydroelastic vibration phenomena are discussed. The sig- 
nificance of the vortex as a forcing function is emphasized. 


The influence of plate vibration upon the early wake is 
essentially eliminated in the investigation by rigidly 
clamping the test plates. 


A prediction of plate profile drag coefficient and of 
Gimensionless shedding frequency, the Strouhal number, is 
made based upon von Karman's vortex street analysis. 


The longitudinal and transverse distribution of local 
mean velocity, turbulence intensity, and static pressure 
were measured in the early wake of three selected test 
plates of varying trailing edge configuration. Turbulence 
energy density spectra were obtained at selected positions 
in the early wake. Profile drag coefficients were computed 
and vortex shedding frequencies measured, and are compared 
with predicted values. 
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I. INTRODUCTION 


A. General Statement of the Problem. 

One of the most vital problems in ship design today, and 
particularly in warship design, is the reduction of noise and 
vibration. Of the many and varied sources of noise and vibra- 
tion in ship structures, one of the more significant is the 
hydroelastic vibration of propeller blades and control sur- 
faces. The phenomenon of "singing" propellers has interested 
workers in several fields since the 1930's, when higher pro- 
peller speeds and the use of airfoil blade sections instead 
of ogival sections became prevalent! 1], Fundamentally, blade 
vibration of some sort has been held responsible in the 
literature for this phenomenon, but writers have differed in 
their analyses of the vibration and in their several explana- 
tions for the cause of this vibration!?-9], 

While propeller in-flow variations, blade back pressure 
variations, and wheel unbalance may indeed contribute to the 
phenomenon, most writers, however, agree that the probable 
cause of singing appears to be the formation and shedding of 
vortices in resonance with the natural frequency of the blade 
and its entrained water. If the vibrations set up by the pro- 
cess of forming of eddies at the trailing edge of a propeller 
blade are in the audible range of frequency, a singing note 
may be emitteq’?!. It has been found that for a given pro- 
peller system singing occurs in a relatively narrow frequency 
band and over a relatively narrow range of propeller rpm. 


While audible singing is generally not encountered in 





larger foil or plate structures such as rudders, hydrofoils, 
submarine diving planes and other control surfaces, vibration 


[10,11]. It appears 


and flutter are significant problems 
reasonable to consider that vortex induced vibration may be 
at least in part responsible for this. 

It has been found in practice that modifications to the 
geometry of propeller blade sections will alter the character- 
istics of singing or eliminate it entirely. Den Hartogl 12] 
and others discuss the procedure of sharpening the leading 
edge of propeller blades and turbine blade runners to eliminate 
undesirable singing. Propeller sections with sharp leading 
edges have, however, been known to sing violently, and several 
authors!1+4. 91 discuss the practice of blunting the trailing 
edges of blades to alleviate this condition. The periodic 
transverse forces due to alternate shedding of vortices incite 
a resonant mode of vibration in some part of the propeller 
blade, and according to Saunders! 2! the remedy is to change 
the shape of the trailing edge in some manner so that the vor- 
tex formation is eliminated or that the vortices and circu- 
lation are stabilized. One simple method is to cut off the 
trailing edge square, so.that the separation points are "well 
defined and fixed in position;" another is to plane the 
trailing edge flat at some angle to the base chord. Various 
other trailing edge geometry modifications are discussed in 
the literature. 

Since practical experience shows that modification of 
the trailing edge will alter the response of a blade section, 


and since it is not reasonable that such modifications would 





materially alter the mass or the internal damping of the 
system’ ’], it follows that it is the forcing function, i.e., 
the shedding of vortices, which has in some way been modified. 

A large body of literature exists on the subject of wake 
flows behind various types of bodies. The majority of it, 
however, deals with the wakes behind bodies of revolution and 
bluff bodies, and almost all of it treats the so-called fully 
developed wake well downstream of the body. In order to con- 
duct a study of trailing edge vortex phenomena at their 
origin, it becomes necessary to investigate the “early wake" 
or the region within about ten characteristic dimensions 
downstream of a foil or flat plate. Such a study should pro- 
perly consider the distribution of the energy present in the 
early wake and the variations of velocity, pressure, and turb- 
ulence intensity in this region. Furthermore, the plate bound- 
ary layer itself and the free shear layer where vorticity 
becomes pronounced may be expected to influence these phenomena. 
These considerations have prompted an investigation by the 
authors of the early wake behind three flat plates of varying 
trailing edge configuration at the variable flow water tunnel 
of the Massachusetts Institute of Technology Hydrodynamics 
Laboratory. 

B. Scope of the Investigation. 

As part of a general investigation of wake mechanics 
conducted at the Massachusetts Institute of Technology Hydro- 
dynamics Laboratory and sponsored by the David Taylor Model 
[13] 


Basin, Huval recently conducted an investigation of the 


early wake behind a fixed flat plate having a semi-cylindrical 





leading edge, a chord-to-thickness ratio of 8.0 and a 60° 
included angle trailing edge. This particular plate was one 
of a series studied by Ippen, Toebes, and Eagleson! 14! for 
the purpose of investigating hydroelastic behavior of flat 
plates with various trailing edge configurations. The authors 
chose three other plates from this systematic series for their 
investigation; their dimensions and geometries are listed 
below in Table I: 

TABLE I 


DIMENSIONS AND GEOMETRIES OF TEST PLATES 


PLATE LEADING EDGE THICKNESS CHORD TRAILING EDGE 


A semi-cylindrical 1/y" Tig ae = 0 
B semi-cylindrical ye De —— 
C  gemi-cylindrical 1/" 2" ae 


These three plates correspond respectively to Plates 12, 11, 
and 15 of Reference [1]. 

The experimental program for this investigation provided 
for the measurement of velocity and pressure distributions 
and distribution of turbulence intensity across the wake at 
several longitudinal positions in the early wake of the test 
plates, and for measurement of the turbulence energy spectra 
at several longitudinal positions on the wake centerline and 
in the vicinity of the paths of shed vortices. 

Plates A and B were originally selected from consider- 
ation of vibrational amplitude responses as a function of 
restraint (Figure 29 of 14); these particular plates repre- 


sented respectively the lowest and one of the highest 





response curves presented. During the course of the investi- 
gation, it became desirable to include an intermediate con- 
figuration, and Plate C was incorporated into the investigation. 
As in Huval's investigation, the authors decided to clamp 

the plate rigidly in the water tunnel test section in order 

to be able to study the wake energy with minimum contribution 
from plate vibrations. In order to facilitate future compari- 
son with Huval's results, it was decided to conduct measure- 
ments at a free stream velocity, Up, in the neighborhood of 


10 feet per second. 





II. REVIEW OF PREVIOUS WORK 


A. Vortex Excited Blade Vibration 

Of the large number of explanations postulated for 
singing propeller blades over the past 30 years, the sugges- 
tion that vortex formation and shedding might be the 
responsible mechanism appears to be one of the most plausible. 
Experimental work conducted in the United States and in the 
Soviet Union in recent yearg! 815] support strongly the con- 
tention that singing occurs when the vortex shedding fre- 
quency is in resonance with one of the natural modes of 
vibration of a given structure. 

Considering that at resonance singing frequency coin- 
cides with vortex shedding frequency, Gonawer!?] presented an 
empirical modification of the Strouhal formulation, suggesting 
that the characteristic trailing edge dimension should include 
a portion of the boundary layer thickness. Introducing this 


factor, he obtained experimentally from his work with vanes 


Spt = u 0.185 (1) 
where b'= ¢ + Q,022ix (2) 


ite 
and where U, is the component of the free stream velocity 
normal to the trailing edge of the body. Krivtsov also 
illustrates!15] the essentially constant behavior of Strouhal 
number for a given trailing edge thickness. 

several writers have been concerned with the fact that 


the singing behavior of identical propellers may be quite 





different. However, when one considers the processes in- 
volved in propeller manufacture and be ldndine, it is not at 
all surprising that two propellers of the same class might 
have slightly different natural frequencies, and that the 
frequencies of the individual blades might vary somewhat. 
Wnkeseer''! points out that slight differences in natural 
frequencies between blades are not at all uncommon. This may 
increase the response or decrease the equivalent damping of 
one propeller with respect to another having blades of iden- 
tical natural frequency. If propeller singing is due to 
self-excited vibration, a small difference in response could 
make all the difference between singing and non-singing. 
Propeller singing has been shown to be dependent upon 


[1] 


water temperature and ambient pressure, indicating further 


the flow dependence of the phenomenon. 
B. Vortex Formation and Wake Perlodicity 

Since Eee ee nea dic work in the latter part of the 19th 
century, a large number of investigators have studied the 
vortex shedding of cylindrical and spherical bodies, and 
more recently of bluff bodies. While this work is of con- 
Siderable interest in theory and in practice, and valuable as 
background, it is not immediately applicable to the case of 
thin flat plates and other blade sections, and will be treated 
Only cursorily. 

The Reynolds number dependence of the Strouhal number, 
Suggested by Raleigh'1°! has been clearly demonstrated 


[17] 


recently for circular cylinders by Kovasznay and by 


Roshkol18) | Roshko presents empirical expressions for Sq 





for two Reynolds number regimes of flow: 


0.212 (1 - 21.2/Reqg) lo < Reg < 150 (3) 


Sq 


and 


Sq = 0.212 (1 - 12.7/Reg) 300 < Reg < 5000 (L,) 


For wakes behind various bluff bodies, Ippen et al. re- 
view the work of several investigators 14! and indicate that 
the general range of Strouhal numbers based upon the largest 
characteristic dimension normal to the direction of the flow 
is from about 0.12 to 0.21, i.e., below the value of 0.2 
applicable to circular cylinders. 

Fage and Johansen in experiments on the early wake bee= 
hind inclined flat plates showed that Strouhal number is 
dependent upon angle of incidence, a, increasing as a becomes 
sma11!19) , Expressing S in terms of the separation between 
free shear layers, S was found to be approximately constant 
and ~ 0.16. Both in Gongwer's experiments and in other ex- 
periments on oval straes’ (a it was found that the effect of 
small angles of incidence was to alter the amplitude of 
vibration without materially affecting the frequency. 

In none of the theoretical treatments of the problem 
is an explanation of the actual mechanism of vortex formation 


presented. The work of Betchov' 2°] 


suggests that a possible 
approach may lie in the treatment of boundary layer oscii- 
lations. 

Von Karman's classical treatment of periodic wakes has 
crystallized and stimulated significant work in this field. 


His stability criterion for a vortex array in a wake has 


theoretical and practical significance and will subsequently 





be utilized to predict shedding frequencies and to analyze 
free shear layer separation. This concept and Birkhoff's 
modification will be discussed further under "Theoretical 
Considerations." Von Karman's treatment makes no attempt to 
Gescribe the mechanism of vortex formation and is in fact 
independent of the mechanism and of the geometry of the 
structure behind which it applies. 
C. Wake Structure 

The similarity principles of von Karman applied to the 
Reynolds equations of motion have proven useful in treating 
axially symmetric or two-dimensional turbulent wakes. Un- 
fortunately, similarity hypotheses apply to a state of 
statistical equilibrium of flow processes relatively far 
Gownstream in a wake, and cannot be considered applicable to 


the early wake. 


[18] [21] 


Roshko's experiments and those of Cooper and Lutzky 
present results of velocity and turbulence distributions in 
the relatively early wake behind cylinders and bluff bodies 
respectively. 

Birknorr! 22] and others point out that contrary to the 
predictions of the Kirchoff wake model, actual large wake 
underpressures exist. 

D. Trailing Edge Geometry 

Ippen et al. point out the effect of trailing edge 
geometry on drag coefficients for various shapes. Schliche- 
ting!23] shows that the drag coefficient can be related to 


the vortex street through momentum analysis, but notes that 


this procedure is applicable only relatively far downstream 





10 
where the static pressure at the measuring station is equal 


to the free stream pressure. The correction methods of Betz 
and of Jones! 241 may be employed in the early wake region, 
and have yielded valid results behind airfoil sections as 
close to the trailing edge as 0.05 chord lengths. 
E. Turbulence 

The complex nature of wake turbulence and turbulent 
(or even periodic) wake energy precludes any detailed dis- 
cussion here. Reference [28] in particular and portions of 
references [29] through [31] are of interest for this inves- 
tigation as are the results of Huval's experiments! 13), 
Ambient turbulence background measurements were made 


[32] and by 


for the water tunnel test section by Perkins 
Huval [13] indicating relative mean turbulence levels of 

about 0.03 and 0.065 respectively. Assuming isotropic tur- 
bulence with a general power decay law, and ergodic process, 
Huval's measurements indicated that the free stream turbulence 
was largely in its initial stage. The character of the free 
stream turbulence spectrum obtained by Huval indicates 
anisotropy in the low frequency end; the high frequency range, 
however, indicates isotropic statistical equilibrium. Huval 
notes that the high frequency spectrum appears to follow von 
Karman's suggested energy variation with ee [13], This 

is in agreement with the Kolmogoroff spectrum law cited by 
Hinzel 28] | 

Hinze points out that in all cases where the mean vel- 


ocity shows a gradient turbulence will be anisotropic. “t+ 


Should then be expected that in regions of a wake where such 





gradients exist anisotropic turbulence will obtain. This con- 
dition is expected to exist in the early wake behind an object. 


[18] notes that the velocity fluctuation in the 


Roshko 
wake of a vortex-shedding cylinder displays a predominant 
frequency which is the shedding frequency; he points out 
further that this predominant frequency will have random 
irregularities imposed upon it. 

Distributions of turbulence intensity, mean velocity, 
and pressure in a turbulent wake are discussed by Hinze. 
Theoretical treatments of these distributions are in general 
predicated upon similarity assumptions which are well known. 
such assumptions, however, are not expected to be justified 
in the early wake behind a structure. Cooper and Lutzky! 211 
found that certain predictions made on the basis of similarity 
hypotheses were not upheld by experimental evidence ina 


region relatively close behind bluff bodies. They note that 


this discrepancy was not resolved. 





III. THEORETICAL CONSIDERATIONS 


A. Analysis of the Vortex Street. 


From consideration of complex potential theory, it can 
be shown! 29) that the complex potential function due to an 
infinite row of equally spaced vortices is given by 


k 
wo =< ln sin = (5) 


from which it can further be shown that for two such rows 


parallel to each other and separated by a distance b, 


sé lin sin : (z - >) - ln sin x (z - $+) (6) 


w= 
where; 
a = longitudinal vortex spacing 
b = transverse vortex spacing 
k = vortex strength ; 
Since u - iv = at , differentiation and separation of 


components yields an expression for vortex velocity in the 
longitudinal direction: 
Uy = of tanh ut (7) 

From stability considerations, von Karman showeal 25:26] 
that the street can be stable only when vortices in one row 
are opposite the interval center between vortices in the other 
row, and then only when the ratio b/a = 0.281. 

Rosenhead (cited in 25) shows that rigid boundaries 
equidistant from the mean centerline of the vortex street 
may affect the stability ratio. If h = the distance between 
the rigid boundaries, then as the ratio a/h varies from 0 to 
0.815, the stability ratio has a unique but continuously 


changing value, going from 0.281 to 0.256. For values of a/h> 





13 
0.815 stability exists for a certain range of values of b/a, 


and for values of a/h > 1.419, the configuration is stable 
for all values of b/a. Birkhoff (cited in 1) predicts a 
stability ratio of b/a = 0.35 from inertia considerations. 

The test section of the water tunnel used in this present 
investigation is rectangular, measuring 9 inches high by 7-1/2 
inches wide. The test plate is mounted vertically in the 
section. The dimension h between rigid boundaries in this 
case is then 7-1/2". Assuming that at the trailing edge 
b = plate thickness, t, and proceeding from the von Karman 
criterion a/b = 0.281 as a starting point, the ratio a/h = 
O.1185. Following a reiterative procedure and assuming 
linear decrease of stability ratio with increase in a/h, the 
Rosenhead treatment indicates that the stability ratio for 
vortex streets in the current investigation would be b/a = 
O.271. In actuality the transverse spacing is expected to 
be greater than the plate thickness due to wake spreading. 
The assumption of linear decrease in b/a is considered 
adequate in the region of a/h concerned. It appears, then, 
that a stability ratio for b/a may lie in the region from 
about 0.27 to about 0.28 for this investigation. 

After the manner of Ippen! 14] and Birkhorr! 27) it can 
be shown that considerations of vorticity diffusion from the 
boundary layer can lead to an expression for the strength, k, 
of the vortices in the wake. If K, be the rate of shedding 
vorticity from each side of an obstacle, then as an approxe 
imation: 


is 2 
Ki = Jcuay = | = uydy = ay (8) 





-— ee ae 


Where uj, is the velocity outside the boundary layer near the 
separation point. If B = K/K,, where K is von Karman's rate 
of vorticity discharge from each side of an obstacle, then 8 
is a measure of the anihilation of vorticity at the trailing 
‘edge. Both Prandtl and Roshko(cited in 1h and 27) have found 
experimentally that B = 0.5 for periodic wakes. Applying this 
value and assuming that uz, = U, we obtain: 

K = 0.25 U4 (9) 
Introducing this into von Karman's expression for vortex 
strength, we obtain: 


ee, WOE i Oma. 
Apa, * Stee (10) 


Introducing this into Eq. (7) we obtain: 


u 
ge (1 = gy) = 0.125 tanh n/a (11) 
from which: 


7 = O05(1%5/1)= 0.5 tannenb/a ) (2925) 


Taking as the stability ratio b/a = 0.271 and operating 
on Eqs. (10) and (11) we obtain: 
uy/U = 0.0955 (13) 
k = 0.277 aU 
From von Karman's relationship between profile drag, Dp» 
and wake vorticity, where Dp = kp> (U - 2uy) ~ gre and from 
Eqs. (13) we obtain: 
Dp = 0.0728 pau? (1),) 


from which: 


Cp, = gpbyr = 0-1456 a/t (15) 





and reintroducing the stability ratio b/a = 0.271: 
Cp, = 0.537 b/t | (16) 


It should be pointed out that Ippen et al. assumed a 
value of K/U2 ~ O.4. which is in consonance with empirical 
results cited by Birkhorrt 27), If, however, the experimental 
result of B = 0.5 is accepted and Eq. (8) is considered valid, 
this would lead to a value of u,/U = 1.265. 

Birkhorf!27] points out that it is a corollary con- 
sequence of the theorem of conservation of momentum that in 
viscous wakes vortex transverse spacing tends to increase; 
this is experimentally observed. For slightly viscous fluids 
it is further noted that transverse vortex spacing is essen- 
tially equal to the wake thickness a little behind the ob- 
stacle. Birkhoff notes that for circular cylinders the 
critical diameter Reynolds number above which the wake first 
becomes periodic is lowered by external turbulence. For flat 
plates parallel to the flow the critical wake-thickness Rey- 
nolds number is about 3.l},/Re] where Re, = eee this is the 
same order of magnitude as in the case of a circular cylinder. 

Following Ippen's procedure, the plate thickness Strouhal 
number for vortex shedding can be estimated from Eqs. (13) 
and (15): 

S_ = AR = (1 ~ a)(t/a) = = 
Introducing Eq. (16) we obtain: 
S, = 0.245 ¢ (17) 





from which: 


SU 
ny = ~~ = OF U (18) 


If wake thickness and transverse vortex separation are taken 
as being approximately equal to plate thickness in the early 
wake, then for t = 0.25" we obtain a predicted vortex shed- 
ding frequency of: 

ny = 11.76 U (cps) (19) 

It should be pointed out that none of the foregoing 

takes into consideration trailing edge geometry. Furthermore, 
these developments do not consider any Reynolds number dee- 
pendence upon either Strouhal number formulation or periodic 


wake structure. 





IV. DESCRIPTION OF EQUIPMENT AND INSTRUMENTATION 


A. The Variable Flow Water Tunnel. 

A detailed description of the variable flow water tunnel 
is given in Reference [33]. Photographs of the tunnel test 
section are presented as Figs. Ia. and Ib.; Fig. Ic.is a 
schematic diagram of the tunnel and its associated flow cir- 
cuit. The hydraulic characteristics of the water tunnel are 
presented in Appendix A of Reference (1h). Tunnel operation 
is possible with any one or a combination of three centri- 
fugal pumps. 

Considerations of minimum mean level of background 
turbulence and consistency of pump operation led the authors 
to conduct all measyrements with Pump III providing the flow. 
Referring to Fig. Ic., the flow was routed from Pump III to 
the right through the 18" distribution circuit, back to the 
left past Pump I (which was isolated from the circuit), and 
up through the water tunnel inlet riser. It had initially 
been planned to operate with both Pumps II and III in oper- 
ation simultaneously, with the stop valve between them shut, 
but flow reversal through Pump II and problems in velocity 
control led to the decision to isolate Pump II from the 
eirguit. 

From the expanding tunnel inlet riser the flow passes 
through four No. 16 mesh screens and a 33" long contraction 
into the test section. The plexiglass test section is 36" 
long and rectangular in section, measuring 9" high and 7-1/2" 


Wide. From the test section, the flow returns to the 








FIGURE Ia. Variable Flow Water Tunnel 
Test Section, Showing Location of Amp- 
lifier and Preamplifier. 





FIGURE Ib. Variable Flow Water Tunnel 
Test Section, Left Side. 
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storage tank. 
B. Test Plates. 

The test plates studied in this investigation are 
described briefly in Table I and in more detail in Reference 
[1h]. They were manufactured from 1/4" aluminum stock and 
carefully machined. The plates spanned the 9" vertical 
dimension of the test section, and presented a 2" chord in 
the direction of the flow. The test plates were anodized to 
prevent corrosion and pitting; this process proved to be 
quite successful. 

C. Test Plate Mounting. 

The test plates were mounted using the spindles illus- 
trated in Figure 6 of Reference [14]. The spindles were 
attached to the plates with set screws such that the center- 
line of the spindles was essentially in line with the 
center of radius of the semi-cylindrical plate leading edge. 
The spindles were clamped to prevent rotation. 

Each of the spindles passed through a roughly rectan- 
gular 1" by 2-1/2" hole in the upper and lower test section 
walls. These holes opened further into larger cavities. The 
upper cavity was originally sealed with a solid cover plate; 
because of the tendency of air bubbles to accumulate in this 
cavity, a vent was subsequently installed in the cover plate. 
The lower spindle passed through a gasketed circular hole in 
the bottom cover plate and was again braced below the test 
section. This arrangement is due to the earlier investigation 


reported in Reference [1]. 
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D, Traversing Mechanism. 


In order to obtain data on transverse and longitudinal 
wake structures with the flow measurement devices to be 
described subsequently, it was necessary to provide for 
traversing in these directions. The traversing mechanism used 
by Huval!13J was employed in this investigation. The in- 
strumentation strut, located just downstream of the test 
plate, passes through the test section side walls through 
slotted blocks which permit a longitudinal traverse. These 
are shown installed in Figs. Ila and IIb. The instrumentation 
strut itself, illustrated in Fig. IIc, can be traversed in 
the transverse direction through the side plates and blocks. 

Longitudinal positioning of the strut is accomplished 
by means of scales graduated in 1/6, inch mounted on the 
slotted blocks. A transverse scale graduated in 0.002 inch 
permits transverse positioning. 

E. Flow Measurement Devices. 

The total and static head tubes and the turbulence 
pressure transducer employed in the wake surveys are described 
in detail in References [13] and [32]. These instruments 
permit measurement of the total and static pressures in the 
wake behind the test plate, and will be treated separately 
below. In addition to these devices, a contraction manometer 
provides for measurement of mean free stream velocity. 

1. Turbulence .Pressure Transducer. 

A cylindrical disc-shaped barium titanate ceramic 
transducer, 0.125" in diameter and 0.100" thick, is mounted 
in a tapered water-filled probe tip illustrated in Fig. III. 





Be 





FIGURE IIa. Left Side View of 
Test Section, Showing Position 
of Test Plate and Transverse 
Strut Positioning Scale. 





FIGURE IIb. Right Side View of 
Test Section, Showing Total and 
Static Head Leads and Transdu- 
cer Cable to Preamplifier. 








FIGURE IIc. Photographs of Traversing System. 
(From Huval, "The Early Wake of a Fixed Flat 
Plate.") 


ey 
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The ceramic is located behind a thin silver-plated brass 
diaphragm. The tapered probe tip ends in a 1/16" diameter 
with a 0.025" hole drilled in it to permit "point" pressure 
fluctuation measurements. Pressure variations are transmitted 
to the ceramic through a water column expanding from 0.025" 
diameter to approximately 0.221" diameter, and through the 
diaphragm. 

Considerable care must be taken to insure that no air 
bubbles are present in the water column. Experience has shown 
that the presence of even minute air bubbles will materially 
affect the response of the system by altering the response 
of the water column. To this end, deaerated water is used for 
filling the probe tip. Reasonable success was experienced 
by pumping deaerated water through the probe tip for several 
minutes with a syringe. 

-t was found earlier that the transducer did not respond 
to static pressure changes 13] Recent improvements, however, 
have improved the response of the transducer and made quasi- 
static calibration possible. This will be discussed further 
in Appendix A. 

2. Transducer Response. 

The instantaneous response of a total head tube aligned 


parallel to the direction of mean flow, in Huval's notation 13] 


ls given by: 


Hie = D +p' + (u+u')e 
¥ MA og 
= p +p' + ue + ouu' +4 ute (20) 
;  & wn 2 


assuming that the terms involving v' and w' do not significantly 





es 
contribute to the total instantaneous head measured in the 
longitudinal direction. 
Now let Cte = the instantaneous output signal, of the 
turbulence gage in millivolts, and 
turbulence gage sensitivity in millivolts 


Rig 


per foot of water. 


Then the turbulence gage output signal becomes: 


ey, = Kygileg = Kt, (2 +p! +u24un' + ul? (21) 
iy Y 2g gz Og 
Now let Gp = preamplifier gain, and 
ey = voltage output of the preamplifier in milli- 
volts, 
Then: 

Cn = GF te (2 + pt + Be +uul + ut?) (22) 

Y 04 Og g 2g 


Furthermore: 


sy Q 
ep = 5) + <5 


a = ae 12 
ey okte (B + us +u 5 
y. 2g 2g 


Referring to the electronics equipment block diagram 
(Fig. V) it is noted that a DC filter (i.e., a 1 MFD capacitor) 
is installed between the preamplifier and the amplifier, since 
only the fluctuating components of the turbulence gage signal 
are of interest at the measuring instruments. Therefore, the 


input signal to the amplifier becomes: 


t 


CD = Cp - en 
= G)Kt, /p' + uu’ +ufe - ) 
Y g 2g (23) 


Letting e, = amplified turbulence signal input to the measur- 
ing instruments, and 


Ga 


amplifier gain, then 
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em = GaGpktg ( 2 +i oe u'? ) (23) 
Y g 2g 

Here em represents the amplified turbulence gage signal in 
millivolts at the input to the measuring instruments (rms 
voltmeters and vibration analyzer). 

The true rms voltmeter output (measured by the thermo- 
couple voltmeter) becomes: 

Vems = Ven? = GaSpktg e tuul eal? aut wi? e (2h) 

Y gz 2g 

Expansion of the terms within parentheses in Eq. (24) results 








in: 
(2 +uu' + u'e@ - A = ueul? + uu'(ul2 - ul@) 
Y 2g ge ge 
+ (u!2 - u'2)2 + Quu'p' 
Age Yg 
+ (u'2 = ut@)p' + p'2 


Yg ye 


The full expression for rms voltmeter output becomes: 

















= GaGrFte Moe Selo + We = ee 
ge ge ge age 
+ ouubeiee ulep" = pe y (25) 
Yg Yg hf 


Now if it may reasonably be assumed that u'<<u, say for 
example that./GTS <0.1 u, and that p'/yY is of the order 
of ule 


pt 


eg 
Significant contribution to the rms signal. Hence: 


ms ~ GaGpki, avy? (26) 


only the first term within the parentheses makes a 


V 
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3, Total Head and Static Pressure Tubes. 

The total head tube employed is in fact a standard Pitot- 
static tube, but only the dynamic tap is utilized. Because of 
the longitudinal pressure gradients present in the wake of the 
test plate and their effect upon measurements, the static taps 
of the total head tube could not be used, and a separate static 
head tube was employed. The pressure taps of the static head 
tube were in line with the tips of the total head tube and the 
transducer probe. See Fig. If. Huva2 L'3) gives a detailed 
description of these instruments. Both pressure tubes are 
connected to mercury manometers. 


4. Contraction Manometer. 


Connected to the upstream and downstream ends of the water 
tunnel contraction are a mercury manometer and a manometer 
filled with blue-colored Meriam No. D-8325 indicating fluid 
(specific gravity = 1.75). These manometers and their piezo=- 
meter taps are described in Appendix A of Reference [14] . The 
contraction manometer was used to duplicate flow conditions in 
the water tunnel on separate days of operation. Calibration of 
this manometer was necessary for the particular flow conditions 
experienced, and is discussed subsequently. 


FF. Electronics Equipment. 


Electronics equipment utilized in the measurements consis- 
ted of an electrometer preamplifier, a KinTel Model 11i1BF 
amplifier, a General Radio Type 762-B vibration S28 oe a 
Ballantine vacuum tube voltmeter and a Sensitive Research Instr. 
Corp. Model A thermocouple voltmeter, and a Dumont Type 333 
dual=-beam cathode ray oscillograph. A Hewlett-Packard Model 


200C audio oscillator was employed in equipment calibration. 
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FIGURE III. Details of Turbulence Probe and Pressure 
Tubes. (From Huval, "The Early Wake of a Fixed Flat 


Plate.") 
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These components are illustrated in Figs. IVa and IVb. 

(At the end of the investigation a General Radio Model 
1445-A vibration analyzer became available; this instrument 
extended the frequency range of spectrum measurements and 
was used in one location in the wake of Plate C for comparison 
purposes.) 

Tt had originally been intended to record the amplified 
signal from the pressure transducer on magnetic instrumentation 
tape for remote spectrum analysis on another wave analyzer. 
Recording difficulties made this impossible during the course 
of this investigation. 

A block diagram of the electronics components is pre= 
sented as Fig. V. 


1, Electrometer Preamplifier. 


The primary purpose of the electrometer preamplifier is 
to provide high input impedance to the first stage of ampli- 
fication. Huva UJ describes the circuitry of this component, 
and points out that its addition to the system should permit 
response to frequencies below i cps. The measured gain of this 
component (see Appendix A) is 0.4 over the majority of the 
operating frequency range. 

2. Amplifier 

The KinTel amplifier used in the system has a very flat 
response over the operating range of frequencies. The design 
gain accuracy with a fixed gain setting is + 1% from DC to 2 ECS. 
The measured response characteristics (see Appendix A) over the 
operating frequency range point this out further. A DC filter 
1s connected to the input of the amplifier. 








FIGURE IVa. Electronics Equip- 
ment: Analysis Components. 





FIGURE IVb. Electronics Equip- 
ment: Preamplifier and Ampli- 
fier Installation. 
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FIGURE V Electronics Equipment Component Block Diagram 
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3. Wave Analyzer. 

The General Radio Type 762-B vibration analyzer is a 
constant percentage band width tunable voltmeter with a 
logarithmic scale. ‘Its operating range is from 2.5 to 750 
cps. For the sharp selectivity position (in which all measure- 
ments were made) the relative attenuation is 3 db at a frequency 
differing by 1% from tuned frequency. Calibration was per- 
formed using a signal from a Hewlett-Packard audio oscilla- 
tor, and is discussed further in Appendix A. A damping 
capacitance was inserted in parallel with the analyzer meter 
to reduce needle fluctuations. 

The Model 1554-A analyzer used briefly at the end of the 
investigation is a similar instrument, but has an extended 
high limit of 25,000 cps for its operating frequency range. 

4, Voltmeters. 





The Ballantine VIVM is a true rms-indicating meter only 
for sinusoidal waveforms. The thermocouple voltmeter, on the 
other hand, will indicate true rms voltages for any signal 
input. 

5. Oscilloscope. 

The Dumont dual=-beam oscillograph was used primarily for 
visual monitoring of the amplified turbulence signal. In this 
respect it proved very useful in determining the presence of 
air bubbles in the turbulence probe tip. Such bubbles attenuate 
severely the higher frequencies of pressure fluctuation. See 


Pigs. Vim, Vib, and Vic. 
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FIGURE VIa. Turbulence Signal FIGURE VIb. Turbulence Signal 
Presentation in Wake of Plate Presentation on Wake Center- 
C, y/t=0.5, Illustrating Pre- line, Illustrating Predomin- 
dominant Vortex Frequency. ant Double Shedding Frequency. 





FIGURE Vic. Turbulence Signal 
Presentation with Air Bubble 

in Transducer Probe Tip, Same 
Location as Figure VIb above. 
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V. CONDUCT OF THE INVESTIGATION 


A. Ambient Turbulence. 

A study was made of the ambient mean level of turbulence 
as indicated by the amplified rms voltage from the pressure 
transducer for various modes of pump operation. As noted 
earlier, Pump III and the previously described flow circuit 
were used in the final investigation; this combination 
yielded one of the lowest mean levels of turbulence. 

It had been decided to operate the tunnel at a flow 
velocity of approximately 10 fps to facilitate correlation 
with Huval's results!13], It subsequently developed after 
final calibration of the contraction manometer that measure- 
ments had been taken at a value of free stream velocity in 
the test section, U,, of about 9.5 fps. 

Energy spectrum measurements were made for the empty 
test section at one longitudinal location on the tunnel 
centerline. 

B. Wake Survey Measurements. 

Pressure, mean local velocity, and turbulence intensity 
distributions across the wake were desired at representative 
locations in the wakes of the three test plates. A coordinate 
System was defined as follows: the x-axis was taken as the 
flow centerline, which was assumed to coincide with the plate 
and tunnel centerlines. The +x-direction was taken in the 
downstream direction, measured from an origin located in the 
transverse plane of the plate trailing edge. The origin of 
the y-axis was taken at the particular longitudinal station 





oD 
where measurements were desired (x/t = 1.50, 3.25, 6.50, and 
9.75). The ty-direction was taken in the transverse (hori- 
zontal) direction toward the right side of the test section 
(see Fig. Ia.). It was decided to measure total and static 
head and turbulence rms voltage across the wake in increments 
of 0.1 plate thickness. Full wake surveys were made from 
y/t = -1.5 to y/t = +1.5 in this manner at longitudinal 
Stations x/t = 3.25 and 9.75. Half-wake Surveys (from the 
flow centerline to y/t = +1.5) were made at x/t = 1.50 ana 
6.50. 
C. Wake Energy Spectrum Measurements. 

Turbulence rms voltage and vibration analyzer measure- 
ments were taken at y/t = 0 (flow centerline) ana y/t = +0.5 
(in line with the plate edge) at each of the longitudinal 
measuring stations mentioned above, for each of the three 
test plates. The empty test section background spectrum 
mentioned in Section V.A. above wag taken at y/t = O and 
x/t = 3.25, 

D. Longitudinal Variation Measurements. 

In addition, it was desired to study the longitudinal 
variation of pressure, velocity, and turbulence intensity 
along the centerline and at y/t = +1.5. To this eng the 
measurements described in Section V.B. were also conducted 
at these transverse locations at the following longitudinal 
locations: x/t = 0.5, 0.75, 1.0, 1.5, 2.0, 2075, 3.25, h.o, 
Paes 6.5, 7.5, 8.5, and 9.75. 


E. Calibration of Equipment . 


Calibration of equipment is discussed in Appendix A, 
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F. Flow Duplication. 

With the foot valve at the discharge end of the water 
tunnel set and left in one position, it was found possible 
to adjust the pump discharge valve and the distribution 
circuit stop valve to duplicate very closely both velocity 
and pressure conditions within the test section. The Meriam 
fluid contraction manometer was used for this purpose. This 
procedure was followed in establishing flow during the 
measurements made for Plates B and C; minor changes were 
made in the foot valve setting during the Plate A measure- 


ments. 
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VI. REDUCTION OF DATA 


A. Mean Flow Velocity, Uo. 

The contraction manometer was used to measure free stream 
velocity. The Meriam blue fluid manometer lent itself to 
more precise reading and was thus used primarily rather than 
the mercury manometer. The contraction manometer was cali- 
brated against a Prandtl-pitot tube inserted through the upper 
tunnel wall at the upstream end of the test section. Velocity 
of the free stream was obtained from the expression: 

Uo = 2.06,y/ Ahy (27) 

where Ahy = Meriam blue fluid manometer difference 
reading in inches. 
B. Local Mean Velocity, u. 

Local mean velocity was computed from the total and 


static head tube manometer readings using the expression: 


5 = Ym 
/ 22 Yw (Aht - Ahs) ape AA (28) 


where: Aht = total head manometer difference in feet 
Ahgs = static head manometer difference in feet 


Az = difference between total and static pressure- 
side manometer readings, in feet 


Because of small variations in values of U, for dif- 
ferent sets of data, it was decided to normalize u with respect 
to the value of u at y/t = +1.5 in every case. 

C. Turbulence Intensity, /ut2 

Following the development set forth in Section IV, the 

turbulence intensity, 7 aoe was calculated from the rms 


turbulence voltage measurements using the expression: 
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Vie. ie (29) 
u GaGpEte 





Turbulence intensity was normalized with respect to the 
value of u obtaining at the point of observation. 


D. Static Pressure, De 





Static pressure was calculated from the static head 
tube manometer reading for each measurement location. It 


was normalized in terms of a pressure coefficient Pp - Do 
$0Uge 
where py, is the free stream pressure and U, is the free 


stream velocity. 


E. Normalized Turbulence Energy Density, E(n)/ u'@. 





Normalized turbulence energy distribution density was 
calculated directly from the thermocouple voltmeter reading 
and the reading of the vibration analyzer at each observation 


frequency from the expression: 


Ein) = _1__ | Rdé(n) (30) 
ule Vragaccy n 


This normalization process eliminates variations in 
absolute magnitude of energy density due to variations in 
u existing over the entire investigation. 


F, Profile Drag Coefficient, Dp. 





The profile drag coefficient was calculated at longi- 
tudinal stations x/t = 3.25 and x/t = 9.75, using the 


Cp, = | |2u U P= Po y 


where the integration was performed numerically by the 


expression: 


trapezoidal rule between the limits y/t = -1.5 and 
y/t = 4#1.5. 
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Applicability of this expression to the early wake was 


checked using the Jones Methoal 233] where: 


1), EE  /EBB)a(¥) 


total pressure at the measuring station 





where: &9 


P> static pressure at the measuring station 
a static free stream pressure 
= total pressure minus static pressure in 
the free stream. 
Integration was performed numerically as above. 
It should be noted that the drag coefficient in terms 


of wetted surface is given by: 


“Dp eG CDi 
for a chord-to-thickness ratio of 16 as is the case for 


these test plates!13), 
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VII. RESULTS OF THE INVESTIGATION 


A. Wake Structure. 

The results of the study of the wake structure of the 
three test plates are presented graphically as Appendix B 
to this report. The salient points evidenced by these plots 
Will be discussed in this section. 

Figs. Bel through B-3 are summary plots, showing for 
each of the three test plates all of the quantities investi- 
gated in the wake surveys. Necessarily drawn to a reduced 
scale, these plots are intended merely to indicate in one 
location the general pattern of the normalized mean velocity 
profiles and longitudinal distributions, the transverse and 
longitudinal distributions of normalized turbulence intensity, 
and the transverse and longitudinal distributions of normal- 
ized static pressure. 

Figs. B-l, through B-6 show for each individual plate 
the development of the dimensionless velocity profile. 

Values of u are normalized in each case with the faired 

value of u obtaining at station y/t = +1.5. Profiles are 
shown for positive values of y/t only. (This is the case for 
all remaining wake plots with the exception of the wake 
symmetry plots.) None of the plates show any distinct in- 
Gication of rapid wake spreading. If the limits of the wake 
are arbitrarily supposed to be the value of y/t for which 
Sot = QO, it 1s seen that for Plate A the limits of the 
wake gradually increase from about y/t = 1.0 to y/t = 1.2 
between the four longitudinal stations. For Plate B the wake 
spreads from about y/t = 0.8 to y/t = 1.23; while for Plate C 
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the wake spreads from about y/t = 1.0 to y/t lil, athe 


velocity defects for Plates A and C at x/t = 3.25 are seen 

to be markedly greater than for Plate C, and the decay of the 
velocity defect is more rapid for these two plates. The 
transverse distributions of velocity defect are essentially 
Gaussian. 

The increased velocity defect at the early stations for 
Plates A and C is illustrated further in Figs. B-7 through 
B-10, which present a comparison of the dimensionless velocity 
profiles at each of the longitudinal measuring stations. 
Close to the trailing edge both of these plates, which have 
relatively sharp separation points, have large defects, al- 
though at x/t = 9.75 all three profiles are essentially the 
same. Plate A, with its 90° reentrant angle trailing edge, 
is seen to be intermediate between the other two plates. 

Figs. B-1l through B-1l present the transverse dis- 
tribution of the normalized turbulence intensity for all 
three plates at each of the longitudinal stations. The wake 
behind Plate C is seen to be the most turbulent in the 
centerline area; that of Plate B is the least turbulent. 
Again Plate A is intermediate. Turbulence intensity decays 
rapidly downstream and at x/t = 9.75 the intensity levels 
are essentially equal and approaching a constant value across 
the wake, indicating a decay to isotropic turbulence. Ture 
bulence intensities appear to approach the free stream level 
by y/t = 1.5 at each station for each plate, and the free 
stream level of turbulence intensity is seen gradually to 


increase as energy is diffused from the wake into the stream. 
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Fig. B-15 through B-18 illustrate the transverse dis- 


tributions of static pressure coefficient in the wake of the 
three plates at each of the longitudinal stations. Plates A 
and C with their more sharply defined separation points are 
seen to have extensive underpressure regions at x/t = 1.53 
Plate B with its rounded trailing edge does not. As we move 
downstream the Plate B distribution becomes partially 
negative, and all three plots move together. At x/t = 9.75 
(see Fig. B-18) the distribution plot for Plate A is seen to 
become markedly negative again. This particular presentation 
is not substantiated by the investigation of longitudinal 
pressure coefficient distribution illustrated in Figs. B-21 
and B-22, and may be due to data error or computation error. 
However, the possibility that this plot is correct must not 
be entirely discounted in view of the unusual trailing edge 
configuration of this plate. The Plate A plot in Fig. B-15 
shows sharp excursions which may be not entirely due to 
scatter. That the sharp drop in pressure coefficient at 
y/t = 0.5 may not be due to data scatter is suggested by the 
fact that Plate B also indicates a smooth decrease in this 
region. It will be seen later that Plate A exhibits a very 
sharp longitudinal pressure coefficient gradient in this 
region and that reverse flow is also observed. 

Fig. B-1l9 presents longitudinal distributions of normal- 
ized mean local velocity at the centerline and at y/t = 1.5. 
In these plots u is normalized with respect to free stream 
velocity, Ug. Reverse flow was observed for Plate A in the 


region from x/t = 0.75 to about x/t = 1.25, and a sharp 
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longitudinal velocity gradient exists. Reverse flow was also 
observed for Plate C close to the trailing edge. 

Fig. B-20 shows longitudinal distributions of turbulence 
intensity for the three plates at the wake centerline and at 
y/t = 1.5. Again Plate C indicates the highest level of wake 
turbulence and Plate B the lowest; Plate A is intermediate. 

Figs. B-2l1 and Be22 illustrate the longitudinal static 
pressure distributions for the wake centerline and at y/t = 
1.5 respectively. Fig. Be21 clearly shows the large under- 
pressure in the wake close to the trailing edge of Plate A. 
Plate C also demonstrates considerable early underpressure, 
while Plate B underpressure is relatively minor and very 
rapidly becomes an overpressure for a short period, as does 
the plot for Plate A. The centerline pressure for Plate C 
is seen to remain negative throughout the longitudinal range. 

Figs. B-23 through B-25 illustrate the essential 
symmetry of the wakes behind the three plates. Although some 
data scatter is seen, it appears that the exact centerline of 
the flow was slightly on the positive side of the established 
axis, fortuitously insuring that the flow centerline was 
actually observed during the half-wake surveys as well as in 
the full-wake surveys. 

B. Profile Drag Coefficients. 

Profile drag coefficients were calculated for the three 
test plates at longitudinal stations x/t = 3.25 and x/t = 9.75. 
Calculations were initially made using Eq. (31), and then 
using the Jones method as indicated in Eq. (32). Results of 


these calculations are listed as follows in Table II. 
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Profile Drag Coefficients 








PLATE x/t Cp, [Ea. (31)] Cp, | Eq. (32)] 
A 3.25 0.443 0.141 
A 9.75 0.378 0.171 
B 3.25 0.519 0.551 
B 9.75 0.318 Oet95 
C D825 0.528 0.196 
C 9.75 0.393 0.308 


The above values of CD) computed from Eq. (31) are 
corrected for the longitudinal pressure gradient by alge] 
braically subtracting from the momentum term a pressure 
coefficient defect term. This term was obtained by numer- 
ieally integrating the defect of the pressure coefficient 
plot where the defect was taken as zero at y/t = +1.5. The 
results compare reasonably with the value of SD) = 0.537 
predicted by Eq. (16). The results yielded by Eq. (32) are 
considered incorrect. 


C. Vortex Shedding Frequency. 


Vortex shedding frequencies were measured by observing 
the peak responses in the normalized energy density spectra 
for the three test plates (see Appendix C). These shedding 
frequencies were observed most strongly at y/t = 0.5, in line 
with the edge of the test plate. The double shedding frequency 
is more clearly noticeable in the spectra taken at y/t = 0, 
on the centerline. Im some cases the triple shedding frequency 
Ls also noticeable. 


Plate thickness Strouhal numbers were calculated for 
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each of the three plates, based upon averaged flow velocities 
and observed shedding frequencies; the results of these cal- 
culations are presented in Table III. 

In order to compare the results of these calculations, 
it was decided to normalize the Strouhal numbers with respect 
to Plate C. It was reasoned that for a small Reynolds number 
range the Strouhal number might be expected to remain 
essentially a constant. Since. variations in Strouhal 
number must be due to variations in the transverse separation 
between vortices, Plate C was selected as a basis for com- 
parison since the square trailing edge provides distinct 
separation points in comparison with either of the other two 
test plates. It was then assumed for the sake of comparison 
that the transverse separation between vortices is indeed 
equal to plate thickness for Plate C. The Plate C Strouhal 
number was then taken as being the constant value for all 
three plates, neglecting small variations due to slightly 
different plate thickness Reynolds numbers, and this Strouhal 
number was considered a vortex-separation Strouhal number, 

Sp = nb/U. From this, the relative values of vortex separ- 
ation were calculated for Plates A and B. These results are 
also listed in Table III, and their physical significance 

is tllustrated in Fig. VII. 
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TABLE IIL 


Strouhal Numbers and Relative Transverse Vortex Spacing 


PLATE sp = & 
A Gee 
B 0. gal 
€ 0.223 


nb 


n U Sp = 7 b 
109.87 9.66 On273 ©. 2565 
114,..87 9.53 OL 2723 ©2222" 
1O0sn00 9.63 On223 0.250" 





FIGURE VII. Relative Transverse Vortex Spacing 


The qualitative representation of vortex spacing in Fig. 


VII illustrates expected behavior in the wake of Plate B; 


the separation points move toward the centerline of the plate 


for a semi-cylindrical trailing edge. When it is realized 


that vortices do not appear at the separation points as dis- 


crete entities, but are gradually formed from vorticity 
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which diffuses out of the boundary layer, even the behavior 
of the vortex locus of Plate A appears plausible. A "cavity" 
of relatively large under-pressure exists immediately behind 
the reentrant wedge, as was pointed out in Section VITZ.A. and 
illustrated in Appendix B, and the flow would be expected to 
be drawn toward the centerline of the wake. The loci of 
vortices are illustrated as gradually diverging from the 
wake centerline in all three cases in conformance with the 
indication of wake spreading seen in Appendix B. 

Observed shedding frequencies were compared with the 
frequencies predicted from Eq. (19). This comparison is pre- 
sented below in Table IV. Again, exact agreement is not 


realized, but a reasonable approximation is indicated. 


TABLE IV 


Comparison of Actual and 
Predicted Shedding Frequencies 








PLATE “Vactual "Vpredicted 
Eq. (19) 
A 110 I Bals 
B es: 112.0 
C 103 (Sa 
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D. Turbulence Energy. 


The results of the investigation of the turbulence 
energy contained at various stations in the early wake behind 
the three test plates are presented graphically in Appendix C. 

Since minor variations existed in the flow conditions 
obtaining for the individual locations, it was felt that no 
precise comparison could be made of absolute magnitudes of 
energy density. To this end, the spectra were normalized as 
indicated earlier. The data presented on these plots is 
independent of the local’ velocity at the point of measure- 
ment. 

The area under these spectrum curves is dimensionless 
and is a measure of the total energy present at a particular 
location in the wake. Integration of the energy density 
spectrum from zero to infinity should properly result ina 
value of unity; this is not the case for any of the spectra 
obtained. This indicates that a certain percentage of the 
total energy present is either not measured or not accounted 
for in some manner. Possible explanations for this fact are: 

(1) The limits of integration are the finite limits 
of the wave analyzer. Clearly some energy must be present 
below 2.5 cps and above 750 cps. Subsequent measurement of 
one of the spectra with a General Radio Type 1554=-A vibration 
analyzer indicates a considerable amount of Whe eee present 
at about 1400 cps, which is possibly attributable to Helm- 
holtz resonance effects. (See Appendix D.) It is not believed 
that any significant amount of energy is present below the 
low frequency limit of the wave analyzer because of the 


small frequency increment under consideration. 
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Comparison of Wave Analyzer Energy Density Spectra 


PIGURE VIIT 
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(2) The spectral measurements made with the wave anal- 
yzer are normalized with respect to measurements made with a 
true rms voltmeter. It is possible that the thermocouple 
voltmeter, being an all-pass energy measuring device and 
indicating a true picture of the total energy present at a 
given location, measures components which the wave analyzer 
meter does not. The wave analyzer is a tunable voltmeter 
designed to measure Fourier components of vibration, and some 
energy components may not be recognized by it. During the 
course of the energy spectrum measurements, simultaneous 
readings were taken with the Ballantine VIVM. On the assumption 
that the foregoing might be correct, the energy spectra were 
also normalized with respect to the Ballantine voltage, since 
1t was reasoned that this instrument is similar to the meter 
installed in the wave analyzer. The results of these calculations, 
Which are not plotted, still indicate an energy defect, but a 
much smaller and more reasonable one. 

Further investigation will be necessary to explain this 
discrepancy fully before any absolute magnitudes can be 
assigned to energy data. 

Because of the foregoing, no magnitude comparison of 
discrete energy represented by the vortex shedding frequencies 
could be made. Im order to examine the vortex energy as a 
function of trailing edge configuration then, the ratio 
A,/Ag was studied, where A, is the area under a discrete 
shedding frequency (or double shedding frequency) peak and 
A, is the total area under the spectral plot obtained. This 
ratio is then a measure of the percentage of energy present 


in the vortex shedding phenomenon for each of the spectra. 
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The ratio AD/As is greatest at any given location for 
Plate B and least for Plate 6; Plate A is again intermediate. 
This fact was noted for both the fundamental and double 
shedding frequencies. The ratio AD/As shows a linear decay 
downstream for both Plates B and C. The value for Plate A 
increases to a maximum in the vicinity of x/t = 3.25 and 
then exhibits essentially a linear decay with x/t. 

Two independent investigations of the ambient level 
of turbulence in the empty test section (with no test plate 
installed) were conducted. The normalized energy density 
spectra for these investigations are presented in Fig. C-15. 
Both spectra were taken at the tunnel location corresponding 
to x/t = 3.25 and y/t = 0. For each investigation the value 
of u measured by the total and static head tubes was found to 
be 9.52 fps. Turbulence intensities for these measurements 
were 0.0352 and 0.0366, considerably lower than the value 
of 0.068 found by Huval when operating at u = 10.2 fps. 

Also presented in Fig. C-15 are the theoretical 


functions suggested respectively by Dryden and von Karman: 


O(n) = Ay and @(n) = Ao 
1 + By@ne (1 + By @n*)5/6 


While reasonable agreement is obtained with either function 

it appears that the experimental data are better represented 

by the Dryden expression in this case. The individual plate 

wake spectra are also well represented by this expression in 

the high frequency region. This is an indication of an 

approach to isotropic turbulence in this region, while anisotropy 


ls indicated in the low frequency plate spectra regions. 
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Figs. C-1 through C-3 illustrate the downstream develop- 
ment of the turbulence spectra for the individual plates 
Along the centerline. Filge. C-4 through C-6 present the 
longitudinal development for the individual plate spectra 
in line with the edge of the plate. Figs. C-7 through C-10 
compare the centerline spectra for the three test plates 
at each of the four longitudinal stations, and Figs. C-11 
through C-14 make a similar comparison for the spectra 
obtained at y/t = 0.5. Inspection of these Figures shows 
the following trends which are worthy of note: 

Between about 400 cps and 750 cps the slope of all 
the spectra is essentially -2 on the log-log plot. At 
frequencies between the double vortex shedding frequency 
and about 400 cps the slope is well represented by -5/3. 

Discrete energy density spikes associated with vortex 
shedding are very pronounced at x/t = 1.50 and at x/t = 
3.25 for both y/t = 0 and y/t = 0.5. Spikes associated 
with the double shedding frequency are pronounced only at 
the centerline. 

Considerable scatter exists in the spectra in the frequency 
range from 2.5 cps to about 20 cps. 

The discrete frequency energy appears to diffuse initially 
to the random components immediately adjacent in frequency. 
These components in turn transfer part of their energy to 
their neighboring components. Thus, as the flow moves down- 
stream the discrete energy spikes are rapidly reduced by 
diffusion of their energy to the random components in a sort 
of "relaxation" process. 


As the flow moves downstream from x/t = 1.50, the 
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turbulence spectra indicate transition from very pronounced 
anisotropic turbulence superposed upon strong discrete 
components of vorticity toward a more isotropic condition 

of statistical equilibrium. Apparently the high frequency 
components achieve equilibrium early in the wake, since the 
portions of the spectra beyond 500 cps do not appear to 
undergo much change as the flow moves downstream. Hinze 
suggests a plausible reason for this early achievement of 
stability in the high frequency end of the spectrum. He 
points out that the flow within the smallest eddies is no 
longer turbulent, but viscous, and that the lower limit on 
turbulent eddy size is determined by dominance of molecular 
effects (28) | In general, the severity of velocity grad= 
lents increases with decreasing eddy size. The rate of 
viscous dissipation of energy associated with large velocity 
gradients then sets the lower limit on eddy size in accordance 
with the principles of conservation of energy. 

Numerical integration of the first four energy density 
spectra taken for Plate C resulted in very large energy defects 
compared to all the other spectra taken. Applying the 
correction suggested above as a possible means of adjusting 
for differences in measuring instruments (i.e., normalizing 
with respect to the Ballantine voltage recorded instead of 
the voltage reading of the thermocouple voltmeter) only 
.Slightly improved the situation. For these spectra A, still 
remained less than 50% . 

Possible explanations for this large discrepancy are: 

(1) Incorrect wave analyzer calibration constant. This 


is discussed further in Appendix A. 
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(2) Abnormally large amounts of energy associated with 
frequencies above 750 cps. As shown in Fig. VIIT, an extended 
spectrum analysis taken with the Type 1554-A analyzer indicates 
about 10% of the spectral energy above 750 cps. However, the 
energy above 750 cps would have to be considerably larger than 
this to account for the discrepancy. 

(3) Presence of an air bubble in the probe tip. ff 
the natural response of the gage-water-column-probe-tip 
system is located at a frequency just above the range of the 
wave analyzer, a large amount of energy will exist at this 
point and will not be measured by the analyzer. ‘It will, how- 
ever, be measured by the rms voltmeter. This is considered 
unlikely, however, since two separate probe tip filling 
operations were involved. 

Fig. VIII indicates that the relative magnitude of the 
four Plate C spectra in question are probably correct, since 
the pattern is essentially duplicated by the spectrum obtained 
with the Type 1554-A analyzer. 

In order to determine whether the probe tip itself was 
subject to vibrations in the highly fluctuating flow field, 
the probe tip opening was plugged with candle wax. fhe 
turbulence probe was then located at x/t = 3.25 and y/t = 0.5 
in the wake of Plate ©. With a flow of U, = 9.55 fps, a 
turbulence intensity of 0.0019 was obtained. With the FKinTel 
amplifier gain increased to 100, the entire spectrum was scanned; 
response was negligible at all frequencies. With the flow secured 
a voltage was still measured with the thermocouple voltmever; 
this background signal, referred to a flow rate of 9.55 fps, 


corresponded to an intensity of turbulence of 0.0006. 
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A graphical presentation of the decay of discrete vortex 


energy is given below in Fig. IX: 
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FIGURE IX Decay of Vortex Energy Downstream. 
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VITZI, CONCIUSIONS 


A study of the graphical presentation of results included 
as Appendices B and C indicates the following trends for these 
test plates at a plate-thickness Reynolds number of 1.9 x 10°: 

(1) Sharp geometry transitions (well defined flow 
separation points) in the trailing edge create a higher level 
of turbulence intensity in the early wake than do gradual 
transitions. Turbulence intensity diffuses rapidly out of 
the wake at a rate proportional to the initial magnitude of 
the intensity and gradually approaches a decay law propor- 
tional approximately to (x/t) 72. 

(2) Strong wake underpressures exist in the early wake 
behind sharp geometry transition trailing edges. 

(3) The energy density spectrum in the early wake 
is anisotropic, particularly close to the trailing edge where 
longitudinal and transverse pressure gradients are greatest, 
and particularly for low frequencies. 

(4) The high frequency end of the energy density 
spectrum in the early wake is seen to approach isotropic 
conditions. This is also true of the ambient tunnel spectrun. 
E(n) is seen to vary essentially as n7* in this region. 

(5) The discrete energy represented by the vortex 
shedding frequency is seen to decay downstream in a manner 
Slmilar to the decay of turbulence intensity. The rounded 
trailing edge geometry favors the formation of vortices of 
higher energy content. 

(6) The magnitudes of the turbulence intensity found 
in the early wake close to the trailing edge are in excess 


of the order of magnitude originally assumed in the 
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development of the expression for turbulence gage response. 
These magnitudes must, therefore, be viewed in a qualitative 
sense only. 

(7) <A significant percentage (as much as 50%) of the 
total kinetic energy of a mixed periodic-turbulent wake is 
associated with discrete vortex shedding frequencies. 

(8) Energy associated with vortex shedding frequencies 
increases initially in the case of the plate with a 90° 
reentrant angle trailing edge (Plate A), and then decays 
approximately linearly downstrean. uve)? found similar 
behavior in the early wake behind a 60° included angle trail- 
ing edge geometry. 

(9) Taking the amount of energy associated with the 
discrete shedding frequencies in the immediate vicinity of 
the trailing edge as a measure of the energy available to 
excite plate vibrations, extrapolation of the trends indicated 
in Fig. IX indicate that Plate B would tend to vibrate most 
and Plate A would vibrate least. This is in agreement with 


results reported by Ippen et ae 


While centerline vortex energy close to the trailing edge 
is greatest for Plate A, wake centerline forces exerted on the 
trailing edge will act equally on either side of the reentrant 


angle; response to such forces would thus inherently be damped. 
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The following recommendations are made for the conduct 
of further investigation of the early wake behind flat plate 
structures: 

(1) The predictability of vortex shedding frequency and 
the behavior of plate-thickness and vortex-spacing Strouhal 
number should be investigated as a function of Reynolds 
number. Vortex spacing should be investigated by traversing 
the early wake with the turbulence transducer and obtaining 
transverse distributions of discrete spectral energy. 

(2) The behavior of profile drag coefficient with var-~ 
lations in trailing edge geometry should be the subject of 
a separate investigation. Techniques and instrumentation 
in this investigation were inadequate for any conclusive 
understanding of such behavior. 

(3) Particular care should be taken to duplicate both 
flow velocity and tunnel test section pressures for absolute 
comparison of successive runs. 

(4) Tunnel wall spindle cavities should be plugged with 
plates drilled to accept plate mounting spindles in order to 
eliminate the possibilities of Helmholtz or stub tube reson= 
ance effects. (See Appendix D.) Any other cavities in the 
tunnel walls in the vicinity of the test plate or transducer 
should be similarly plugged or faired. 

(S) Frequent periodic calibration of measurement equip- 
ment and components, particularly when they are battery 
powered, is mandatory when accurate data is to be obtained 


over a period of time. 
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(6) Accurate analysis of energy density spectra requires 
accurate calibration of analysis equipment. True rms indication 
of both total energy level and energy density as a function 
of frequency is desirable. The same instrument or accurately 
cross-calibrated instruments should be used to measure these 
quantities. 

(7) Consideration should be given to utilizing frequency 
modulated tape recording equipment to permit determination 


of autocorrelation functions for the energy density spectra. 
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CALIBRATION PROCEDURE 
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1. Turbulence Pressure Transducer. 

a. General. 

The turbulence gage used in this investigation was the 
same as that employed by Huva1!13), A quasi-static calibration 
procedure similar to that used by Huval was used. However, 
the procedure was simplified as a result of repairs to the 
transducer made by Hydrodynamics Laboratory personnel which 
resulted in a greatly increased gage voltage decay time 
constant. 

Huval found the turbulence gage sensitivity to be 63 
millivolts per psi, a value substantially lower than the 
theoretically predicted sensitivity of 100 millivolts per 
psi. Furthermore, his calibration! 23] showed the gage time 
constant to be in the vicinity of 0.6 second, considerably 
lower than that to be expected. Subsequent to Huval's in- 
vestigation the sensitivity of the gage decreased markedly. 
Hydrodynamics Laboratory personnel then disassembled the 
gage and found a small amount of foreign matter on the 
ceramic face. The ceramic was carefully cleaned and re- 
assembled. Tests performed prior to the present investi- 
gation disclosed that the sensitivity of the gage was in- 
creased and the gage time constant was greatly improved. 

b. Description of Calibration Procedure. 

Fig. A-l1 shows schematically the turbulence gage cali- 
bration equipment used in this investigation. 
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constriction 


mercury manometer 
solenoid 
valve 





pressure 
chamber 





pressure bulb 






transducer leads 


FIGURE A-1 Gage Calibration Schematic 


The turbulence gage with probe tip installed was in- 
serted in the tube connected to the side of the cylindrical 
pressure tank as shown. A differential pressure was built 
up in the mercury manometer by means of the hand-operated 
pressure bulb. The clamp in the manometer pressure tube was 
set to give sufficient damping constriction to reduce pressure 
oscillations. The solenoid valve at the top of the pressure 
chamber was opened by means of a trigger circuit which simule 
Ceneously triggered the horizontal sweep on the oscilloscope. 
An oscilloscope camera was used to record the indicated gage 
response voltage. See Fig. A~2. 

A simpler method of calibration, eliminating the use of 


the oscilloscope camera, was employed for the majority of the 
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turbulence gage sensitivity calibrations in this investigation. 
In this alternative method, the oscilloscope sweep rate was 
increased and the horizontal sweep selector set to "recurring 
sweep." The vertical deflection of the horizontal sweep at 
time of pressure release was measured directly on the face of 
the oscilloscope, manually recording on a strip of paper. 

One series of photographic calibration runs was made immediate- 
ly after a series of runs using the simplified procedure. 
Reproducibility between the two methods was excellent, both 


yielding the same value for gage sensitivity. 


c. Development of the Expression for Turbulence Gage 
Sensitivity and Wine tonstants. 


The following definitions will be used in this develop- 
ment : 
Keo = turbulence gage sensitivity 


Ktp = transducer-preamplifier system sensitivity 


Gp = preamplifier gain 

Ty = voltage time constant 

Tp = pressure time constant 

tye = 1/7, 

oan [ie 

p = instantaneous pressure at transducer 
Po = initial calibration pressure 

V = instantaneous voltage at oscilloscope 


As indicated by Huvai £13) , the sudden release of 
pressure from the tank, assuming a suitably damped system, 


may be represented by the following expressions: 





p= poe -P* , for t >O 


Po , tom 6 <0 (A-1) 


p 


The incremental change in voltage generated by the tur- 
bulence gage is the sum of two components: 


AV] = -Kgpdp (pressure-generated component ) 


-CyVAt (voltage decay component) 


AV> 
Hence: AV = AV] + AVo = -KypAp - CyVAt (A-2) 


Rewriting Eq. (A-2) as: 


(> 7) at = “Kep (42 )At - Cat 


and taking the limit as At— 0, we obtain: 


Gv-_ _ dp _ = 
at Kip at rad (4=3) 
Substituting p= poe -P* into Eq. (A-3) we obtain 
av _ =-Cyot 
aE = KypPoCpe PY - CyV (A-l) 
Defining V, = Ky pPo? we see that: 
av _ -Cyt 
at Vane p -CyV (A-S) 
Multiplying both sides of Eq. (A-S) by ev’: 
acvt av Cyt (Cr=c. )t 
VY ae + CyWe Y = Vic,e.°V-P (A-6) 
Then: 
d C a 
at (Ve vt) = Volpe Cpt (A-7) 
Integrating: 


Cyt _ Vo © 7 
veCve = 2B elCv-Cp)t 4 ¢, (A-8) 


ve*P 
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where Cy is the constant of integration, From the initial 


condition that V = O at t = O, we find that: 
VoC 
c ieee (A-9) 
2 Ct, 


Then: 


Vol = “ 
= gk (oot 4s Cyt (A-10) 
V 


Pp 
In analysis of the oscilloscope calibration photo- 
graphs, it is convenlent to express V in terms of V,, the 


maximum voltage occurring at time t,. (See Fig. A-2). 


(ave = +, = 9 


From: 


we obtain: 





-Cyot Cy _-Cyt 
e “Pm = e “-v"m (A-11) 
Cp 
Also: a 
Vn = VoCp (<-°rtn -e-Cv'n) (A~12) 
Cy~Cy 


Substituting Eq. (A-11) into Eq. (A-12) and simplifying: 


Vm = Voe" vim = Kt pPoe” tiv tv 


(A-13) 
solving for Kip in terms of Vm, we obtain: 
Ktyp = Vin vty (A-1)) 
p Po 


Before going on to the development of expressions for 
the time constants Ty and Tp, it should be noted from the 


foregoing expression that if tm < Ty, then: 


Vin 
K os ait 
eP hg (A-15) 
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This yields a means of obtaining an adequate approx- 


imation for K without accurate knowledge of either Ty 


tp 


‘TE. 
or p 


Substituting Eq. (A-13) into Eq. (A-10), the following 


expression for V may be found: 


tt t t 
V= chan GE) -e ~ Th) (A-16) 


If it can be assumed that Tp << Ty and t,<<Ty, the 
second term within brackets in Eq. (A-16) becomes negligible 
with respect to the first term for t >> tp. 


Then: ae, 
V Vme lv ; for t>> tm (A-17) 
Hence: 
x Cet 
Ty = Tn(Vm/Ty 3 evaluated at t>> Cn (A-18) 


To obtain a means of evaluating Tp, evaluating dV/dt 
from Eq. (A-16) above at t = O results in: 


(at) _ Vm 
7); -9 


and: 
Tp = Tav7atT $ evaluated at t = O (A-19) 


This allows evaluation of Tp from knowledge of the 
initial slope of the V versus t plot. 
dad. Results of Pressure Transducer Calibration. 


A series of seven photographic calibration runs were 





68 
made using various values of initial pressure p, and using an 
oscilloscope horizontal sweep rate, calibrated, of 500 milli- 
seconds per inch. A DuMont Oscillograph-Record camera Type 
297, using an f-stop of f/5.6 with a 5-second exposure, the 
sweep being triggered 2 seconds after the shutter was opened. 
Polapan 3000 Type 17 film was used. The time to reach maximum 
voltage, Cn? was about 30 milliseconds on each run, as com- 
pared with a value of about 50 milliseconds reported by Huval. 


(1) Pressure Decay Time Constant, Tp. 
Eq. (A-19) was used to determine T,.. Values ranged 


p 

from 11 to 18 milliseconds, with five of the seven values 
being in the range of 12.5 to 1) milliseconds. An average 
value of Tp = 13.7 milliseconds was obtained. 

(2) Voltage Decay Time Constant, Ty. 

A plot of 1n(Vm/V) versus (t-t,) was used to 
evaluate Tp. See Eq. (A-18). A range of values from 2.) 
seconds to 11.l, seconds was obtained, with five of the seven 
runs in the range l.l, to 9.0 seconds. 

Much of this scatter was due to inability to obtain an 
absolutely stable reference voltage on the oscilloscope prior 
to triggering the pressure release. Also, extreme sensitivity 
of the gage and preamplifier to small unavoidable disturbances 
during the 2-second period of the horizontal sweep may have 
contributed to the scatter. 

The two extreme values were disregarded, resulting in 


an average from the remaining five runs of T, = 6.) seconds. 


p 
However, even using the smallest value of Ty and the 


largest value of Tp» it is seen that: 
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di: tm 
Or = 0.01 
1 0.0075 and Tr 25 


This shows that the approximation of Eq. (A-15) is 
valid for determining Ktp and justifies using the simplified 
procedure described in Section A.1l.b above. 

(3) Turbulence Gage Sensitivity. 

Turbulence gage sensitivity was determined from 
Eq. (A-15) and the relationship: 

Keg = Kep/Gp 

Seven separate series of calibration runs were performed 
during the course of this investigation between 22 February 
and 11 April 1961. The dates of the runs are presented in 
the legend of Fig. A-=-3. 

During the observations made on 6 April, the KinTel 
amplifier was introduced into the calibration circuit to 
provide higher output voltage levels. Amplifier gain setting 
was 30. It was felt that this procedure would serve as a 
further check on the sensitivity of the entire measuring 
system. Due to the quasi-static nature of the calibration 
procedure, the DC filter between the preamplifier and the 
amplifier was bypassed. The average value of Key determined 
for these runs was about 20°/o lower than the February cali- 
bration. 

As a check on possible time-deterioration of sensitivity 
the February procedure was duplicated on 11 April; these runs 
produced a value of K,, about 4°f lower than the original 
February value. It was concluded that the amplifier must 


have attenuated the quasi-steady gage response. 
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FIGURE A-3 Transducer-Preamplifier sensitivity Calibration 
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The amplifier had already been calibrated using an AC 


signal source and found to have excellent gain response. The 
amplifier was now calibrated using a DC dry cell through a 
potentiometer as a signal source. With gain setting of 30, 
actual gain was measured to be equal to about 27. When the 
ratio of gain factors was applied to the gage calibration 
constants obtained with the amplifier in the circuit, the 
same °/o variation from the initial calibration sensitivity 
was obtained. It is considered probable that this variation 
is due to data scatter. 


Finally, an average value of K,. = 16.0 mv/ft. of water 


p 
= 36.9 mv/psi was obtained. This corresponds to a value of 
Kio = hO mv/ft. of water = 92.) mv/psi, as compared to the 
theoretical prediction of Kee = 100 mv/psi reported by 
Perkins and Eagleson! 32!, 


2. .Wave Analyzer. 

Both wave analyzers used in the investigation were 
calibrated using a "pure tone" sinusoidal input signal from 
the Hewlett-Packard audio oscillator. Rms input voltages 
were measured with the Ballantine VTVM. All calibrations were 
performed with the vibration analyzer selectivity in the 
"sharp" position and with vibration analyzer sensitivity set 
at 10. (All data were taken with these settings also.) 

a. General. 

As indicated by Perkins and Fagleson! 32), the General 
Radio Type 762-B Vibration Analyzer output voltmeter reading 
squared i85 proportional to the total energy level of the 


frequency components of the input signal within the narrow 





band to which the analyzer is tuned. 


b. Wave Analyzer Response. 


No An 


Baw = E(n)dn (A-20) 
nyt 


For an input signal with a continuous spectrum of 
frequencies, if E(n) may be considered constant over the 


narrow bandwidth of the wave analyzer, (A-20) becomes : 
Boy = E(ng)An = Cy E(ng)ng (A-21) 
From the statement made in paragraph (a) above, there 


results: 


Ee = CoRa* (A=22) 


Therefore, combining (A-21) and (A-21) : 


2 

ig) = ee (A-23) 
Cy 

where Cy, = 


ce Wave Analyzer .Calibration Procedure. 


The purpose of the wave analyzer calibration is to dee 
termine the constant C, in (A-23) . 

One means of obtaining this constant is to set the wave 
analyzer at a given frequency ny» then to apply an input 
sinusoidal signal from a sine-wave generator. The input 


Signal frequency is increased in small increments from a 


value (well below n,} where negligible wave analyzer response 


is noted, to a value well above n, where again negligibie 


ie 





{@, 
analyzer response results. The input rms voltage €,, as 
measured by an A.C. rms voltmeter, should be kept constant 
throughout this procedure at a value predetermined to give 
nearly full scale deflection at maximum wave analyzer response. 
ey is recorded and constantly monitored. At each input fre- 
quency n, n and Rd are recorded. Then Cy may be determined 


from the relationship: 
12 8 
1 Ra? 
C, = = | No dn (A=2h) 
a 
iL o 
The foregoing procedure should be repeated for several 
values of n, over the range of frequencies which are of 
interest. By this means, any frequency dependency of Cy may 
be detected. 
An alternative procedure (introducing a small error) 
is to set the input frequency at a constant value nj» then 
shift the wave analyzer frequency in small increments from 
the point of negligible response on the lower side of ny tO 
the point of negligible response on the upper side of n° 
Then C., may be approximated from the relationship: 


a Ra“ 
Cy = 5 | - dn, ({A-25) 


The latter method was used throughout this investigation, 
due to instability of the sine wave generator output voltage 
introduced by shifting frequencies. Numerical integration 


was used to approximate the value of the integral. 
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d. Results of Wave Analyzer Calibrations. 
It was noted during this investigation that the wave 


analyzer calibration constant was dependent upon plate vol- 
tage (supplied by two 5 volt batteries). 

This was true even when a check of plate battery voltages 
showed the plate voltage to be within the satisfactory 
operating range. Therefore, frequent wave analyzer calibration 
is necessary in order to obtain accurate results. Particular- 
ly, the wave analyzer must be recalibrated whenever the 
plate supply batteries are changed. 

It is possible that large discrepancies in the areas of 
four energy spectra for Plate C may be attributed to failure 
to calibrate the wave analyzer immediately after a change of 
plate supply batteries. 

On March 28, 1961, wave analyzer plate supply batteries 
were replaced. The turbulence spectrum measurements were 
commenced and all eight energy spectra for Plate C were 
taken between March 28 and March 31, when the wave analyzer 
was calibrated. The wave analyzer constant was found to be 
in good agreement with the two previous values determined. 
However, it was later discovered that the integrals of the 
energy density spectra were very low for the four spectra 
measured on March 28 and March 29 compared to those taken on 
March 30, and, in fact, compared to all the spectra taken for 
the other two plates. At present, lacking a better reason 
for this discrepancy, it is conceivable that the wave analyzer 


constant was very low immediately after the battery replacement 





> 
for some unexplainable reason. In evaluation of this possi- 


bility, it must be noted that a very low value compared to 
any previously obtained would have been necessary to explain 
the discrepancy, and the constant would have increased with 
analyzer operating time, contrary to the trend previously 
experienced. 

The following values of Cy were obtained during this 
investigation. Only average values are listed, since each 


calibration showed little frequency dependence of the constant. 


ea 
= a “te ae 
3-6-61 to 3-10-61 5,12 x 105 
3=10=61 to 3-28-61 Sorex 165 
3-28-61 £6 3-31+61 5.351 x 105 


3..-Amplifier and Preamplifier. 
The electrometer preamplifier and the KinTel amplifier 


were both calibrated using a sinusoidal input signal from 
a Hewlett-Packaré Model 200C audio oscillator. 

Calibration of these components was conducted before 
and after data taking; no time variation of response was 
observed. The initial calibration was performed for both 
instruments by measuring both input and output signals with 
the Ballantine VIVM. This necessitated removing the volte 
meter from the input circuit to measure output voltages from 
the component being calibrated. While this procedure is not 
considered sound practice, the results of the initial cali- 
bration were corroborated by subsequent calibrations, The 
low frequency limit of the audio oscillator used as a sigral 


generator is 16 cps. Hence, the equipment could not be 





76 
calibrated for frequencies below this limit. The initial 
calibrations were conducted from 16 cps in harmonic increments 
of 60 cps to an upper limit of 720 cps with continuous oscil- 
loscope monitoring. Gain factors were computed for the pre- 
amplifier, and for the amplifier at gain settings of 20, 30, 
50, 70, and 100. 

The final calibrations were conducted using the Ballan- 
tine VIVM for input measurements and the thermocouple volt- 
meter for output measurements. Simultaneous measurements were 
taken without disturbing the calibration circuitry. The two 
meters were previously calibrated against each other with a 
sinusoidal input signal over the entire frequency range. 
Finally, the preamplifier and amplifier were calibrated to- 
gether as a system from 16 eps to 5000 cps with amplifier 
gain settings of 20 and 30 (i.e., the gain settings actually 
used during data taking). 

The results of these calibrations are presented in Fig. 
A-4 as a plot of relative attenuation or loss in gain in db 
versus frequency. 

The measured gain of the preamplifier is seen to exhibit 
some attenuation above 350 cps; the response of the amplifier 
is essentially flat up to 750 cps, which is the upper limit 
of the Type 762-B Vibration Analyzer. The preamplifier gain 
was taken as constant over the operating range of the ine 
vestigation at a value Oi 0.40, which was the arithmetic 
mean value. The amplifier gain was G, = 20 with a gain 


setting of 20 and Gz = 30 with a gain setting of 30. 
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FIGURE A=-4 Component and System Gain Response 





The electrometer preamplifier was calibrated and 
operated throughout with a screen grid bias of +2.0 volts. 
A circuit diagram of the electrometer preamplifier is pre= 


sented in Fig. A-5. 
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FIGURE A-5 Electrometer Preamplifier Circuit Diagram 
(From Huval, "The Early Wake of a Fixed, Flat Plate.") 
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FIGURE B-2 Wake Structure Summary Plot, Plate B 
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FIGURE B-3 Wake Structure Summary Plot, Plate C 
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FIGURE B-4 Transverse Velocity Profile, Plate A 
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FIGURE B-6 Transverse Velocity Profile, Plate © 
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FIGURE B-10 Comparison of Velocity Profiles at x/t = 9.75 
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FIGURE B-11 Turbulence Intensity Distribution at x/t = 1.50 
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FIGURE B-12 Turbulence Intensity Distribution at x/t 
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FIGURE B-15 Turbulence Intensity Distribution at x/t 
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FIGURE B-14 Turbulence Intensity Distribution at x/t = 9 
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FIGURE B-15 Pressure Coefficient Distribution at x/t = 1.50 
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FIGURE B-16 Pressure Coefficient Distribution at x/t = 3.25 
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FIGURE B-17 Pressure Coefficient Distribution at x/t = 6.50 
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FIGURE B-19 Longitudinal Dimensionless Velocity Distribution 
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FIGURE B-20 Longitudinal Turbulence Intensity Distribution 
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FIGURE B=-21 Pressure Coefficient Distribution at y/t = 0 
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FIGURE B-22 Pressure Coefficient Distribution at y/t = 1.5 
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FIGURE B-23 Wake Symmetry Plot, Plate A 
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FIGURE B-24 Wake Symmetry Plot, Plate B 
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FIGURE Be25 Wake Symmetry Plot, Plate ¢ 





—_— lage 


1.0 


0.9 


-0.8 | 


0. 6 : 


{025 


0.4 


0.3 


+y/t at x/t 
~y/t at x/t 
+ty/t at, x/t. 
ay /t at x/t 


pb be O 6 
“NMA 
STN 


Heo 


f..| 


ot 


™m% SC“ 
a4 


Aeoh 


4 Re 


a4 fA 00 


Uu 


FIGURE B-25 Wake Symmetry Plot, Plate c 
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APPENDIX C 
NORMALIZED. ENERGY DENSITY SPECTRA 
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FIGURE C-1 Energy Density Spectra at y/t = 0, Plate A 
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FIGURE C-2 Energy Density Spectra at y/t = 0, Plate B 
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FIGURE C-3 Energy Density Spectra at y/t = 0, Plate © 
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FIGURE 0-4 Energy Density Spectra at y/t = 0.5, Plate A 
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FIGURE 0-5 Energy Density Spectra at y/t = 0.5, Plate B 
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FIGURE C-6 Energy Density Spectra at y/t = 0.5, Plate C 
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FIGURE C-7 Energy Density Spectra at y/t = 0, x/t 
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FIGURE C-8 Energy Density Spectra at y/t = 0, x/t = 3.25 
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FIGURE C-9 Energy Density Spectra at y/t = 0, x/t = 6.50 
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FIGURE C-10 Energy Density Spectra at y/t = 0, x/t = 9.75 
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FIGURE C-11 Energy Density Spectra, y/t = 0.5, x/t = 1.50 
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FIGURE C-12 Energy Density Spectra, y/t = 0.5, x/t 
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FIGURE 0-15 Energy Density Spectra, y/t = 0.5, x/t = 6.50 
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FIGURE 0-14 Energy Density Spectra, y/t = 0.5, x/t = 9.75 
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FIGURE C-15 Energy Density Spectrum, Empty Test Section 
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APPENDIX D 
HELMHOLTZ RESONANCE EFFECTS 
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1. Plate Spindle Cavities. 

During measurements incident to determination of the 
characteristics of the ambient tunnel turbulence spectrum 
a possible Helmholtz resonance effect was observed. With 
the desired flow established, the amplified turbulence in- 
tensity voltage was observed to be O mv rms. Monitoring 
the entire available spectrum with the vibration analyzer 
Showed no significant peaks. The test plate had been re- 
moved from the tunnel test section; no plugs had been placed 
in either of the spindle cavities. 

After several minutes the turbulence voltage level began 
to rise, reaching approximately 90 mv rms, and a discrete 
quantity of energy was observed with the vibration analyzer 
at about lO cps. The remainder of the spectrum indicated an 
essentially constant level, down about 6 db from the peak 
reading. The rms voltage level gradually increased and the 
peak frequency was gradually noted to decrease, reaching as 
low as 33 cps. 

At this point a large accumulation of air bubbles was 
noted in the upper spindle cavity under the top cover plate. 
The vent was opened, and as the bubbles vented off, the rms 
voltage level decreased to about 1.0 mv and the discrete peak 
disappeared. 

This phenomenon indicated the possibility of Helmholtz 
resonance to the authors. Clearly only the air bubbles were 
causing the low frequency energy, but even a water-filled 
cavity might be expected to resonate and contribute energy to 


the flow. Even with the test plate installed and the mounting 
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spindles passing through the upper and lower tunnel walls, 


cavities of complicated geometry exist. See Fig. D-1 below. 





FIGURE Del. Sketch of Spindle Cavity 


The fundamental resonant frequency of a Helmholtz type 


3331 


cavity can be estimate using the electrical analogy of 
a simple tank circuit. For a cavity whose dimensions are 
known: 
= surface area of the cavity opening, 
V = volume of the cavity, 
1, = effective length of the "neck" of the cavity, 
approximated by actual length + 0.8/8, and 
c = sonic velocity in the medium, taken as 1:79) 


fps in fresh water, uncorrected for rigidity 


of cavity "neck." 
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Then the Helmholtz resonant frequency, ny, is approximated by: 


Ny = Ts iv (D-1) 


and in this case from actual measurements: 


na B/E, = 250 om 


It is well known that the theoretically derived 
frequency is an upper limit and that in general any number 
of effects can reduce the derived frequency. Furthermore, 
Morse points eut' 32] that the lumped circuit approximation 
is less valid for acoustical systems than it is for electrical 
circuits because of the shortness of the wave lengths in- 
volved. He notes that Eq. (D-1) is of questionable validity 
for frequencies above 1,000 cps, for when a system is longer 
than a half wave length it begins to Behave like a trans- 
mission line and the wave equations become involved. 

2. The Turbulence Probe Tip as_a Helmholtz Resonator. 

The above discussion led naturally to a consideration of 
the turbulence probe tip as a Helmholtz resonator. From the 
interior dimensions of the probe tip as cited in Reference (13), 
and using Eq. (D-1), it was found that for the turbulence 
probe tip: 

Ny 2 7060 eps 

It is considered probable that signals above this fre- 

quency will be attenuated at the rate of approximately 20 db 


per decade. 
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APPENDIX E 
REPRESENTATIVE SAMPLES OF DATA 





et 


SAMPLE WAKE SURVEY DATA 





Pare Ce. 2 /i =.9.75 DATE: 3-28-61 
AMPLIFIER GAIN: 30 


CONTRACTION MANOMETER: Ahy = 21.30, Ahig =p 


VALVE INDICATORS - FOOT VALVE: 1.9 
PUMP DISCHARGE: 22.5 
CLRCULT SPGP: LARS: 


onp as 2BP 82S 2 2 FP SP He Ss HP SF FSF SFP FF SSF PSP TF ES HP HP SH Fe TT Tse FO el EellUDSllCLDl 


TURBULENCE VOLTAGE (RMS) TOTAT, HEAD STATIC HEAD 
Thermocouple Ballantine Pressure Atmos. Pressure Atmos. 

Side Side Side Side 

y/t mv mv "Hg "Hg "Hg "Hg 
41.5 110 93 20.73 25.86 20.55 24.47 
+1.4 115 102 2Om(4+ ~25.65 20.56 o4 47 
it. 5 120 105 20815 25.84 20.56 24.46 
+1.2 125 110 20.76 25.84 20.56 4.46 
+1.1 135 122 20.78 25.82 2025 (. 24,45 
+1.0 139 127 20.80 CD a COs ( 24.45 
+0 .9 145 130 20.82 ls). 1d 20.57 24.45 
+0.8 151 132 20. 83 25.74 20.57 24.45 
+0.7 153 136 20.84 els (i. 20.58 24,44 
+0.6 156 ise 20.86 25.72 20.58 24, 44 
+0.5 157 138 20 107 25.71 20.58 24,44 
+0.4 155 {42 20.88 ean 20.58 24,44 
+0.3 155 141 20.90 25.68 20.59 24.43 
+0.2 157 ee) 20.90 25.68 20.59 24.43 
+0. 1 154 142 20. 90 25.68 20.59 24.43 
O 150 138 20.90 25.68 20.60 24,42 
-0.1 155 143 20.90 25.68 20.59 e443 
-0.2 157 145 20.88 25.10 20.59 24,43 
-0.3 154 140 20.86 25.73 20.59 24.43 
-0.4 149 142 20.86 25.73 20.59 24.43 
=0.5 149 139 20.82 25.76 20.57 24.45 
-0.6 150 138 20.80 25.78 20.56 O4 46 
-0.7 147 132 20.80 25.78 20.56 24 46 
-0.8 142 130 Bore t 25.81 20.56 O4 47 
“0.9 135 123 20.76 25.83 20.56 O4 47 
-1.0 133 123 ZO > 25.84 20.56 4 AT 
-1.1 122 110 20.74 25.85 CO. o> O4 47 
-1.2 117 103 20.74 25.85 20.55 24 47 
-1.3 115 95 20.73 25.86 20S O4 AT 
-1.4 100 oie 20.72 25,86 20755 O4 47 


=1.5 Tf 82 COR em e tO, 20.56 24.48 








SAMPLE ENERGY SPECTRUM DATA 
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DATE: 3-29-61 


PLATE C 
AMPLIFTER GATN: 20 


AVERAGE RMS TURBULENCE VOLTAGE ~ THERMOCOUPLE: 151 mv 


/ Ga See 





ANALYZER SENSITIVITY: 


10 





ANALYZER SELECTIVITY: SHARP 


BALUANTINE: 


ANALYZER CALIBRATION CONSTANT: 5,31 x 10° volts” 


11 


1 


CONTRACTION MANOMETER: Ahy = 21.35, Anup =O 


Ahtotal = 25.62 = 20.95 


A 


mv 


hetatic = 24.36 - 20.66 


FREQUENCY ANANYZER FREQUENCY ANALYZER FREQUENCY ANALYZER 


(cps) READING 


OO ON AWW UIP 


n 


S1) 


READING 


Rd 


ADO DDO ONM OFUFEAVAHNOWWO FU 


pe ee ee ee ee ee ee ee ee ee ee re le 


20 


(cps) 


n 


70 

( 

80 

85 

90 

=e 

98 
100 
102 
104 
106 
108 
110 
HS 
120 
130 
140 
150 
160 
170 
180 
190 
195 
200 
285 
210 
220 





READING 


Rd 


yey) 
25 
28 
31 
38 
60 
cee 
92 


(cps) 


n 


250 
240 
250 
260 
2/0 
280 
290 
500 
320 
340 
360 
380 
400 
420 
440 
460 
480 
500 
B20 
540 
560 
580 
600 
650 
700 
750 


Rd 


28 
26 
24 
25 
24 
27 
26 
27 
25 
25 
25 
25 
25 
24 
23 
22 
22 
22 
22 
22 
22 
O1 
21 
20 
20 
20 


ee nn a | | | | | | - | ; a. | 


NOTE: Analyzer readings from n 


2.5 through n = 


50 are time 
averages of a number of individual readings taken at fixed 
time increments in an attempt to smooth out data scatter. 
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